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Abstract 
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inclusive semileptonic r-decay and study the renormalization scheme depen- 
dence. We argue that the renormalization scheme ambiguity is considerably 
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I. INTRODUCTION 



The r-lepton is the only known lepton massive enough to decay into hadrons. Its semilep- 
tonic decays are very convenient for studying strong interactions at low energies, because 
nonperturbative contributions are rather small and one can apply the standard methods of 
quantum field theory. The measurement of the quantity R T allows one to extract the value 
of the strong coupling constant at a low energy scale. The comparison of this value with ones 
obtained at higher energies is an important test of the applicability of QCD perturbation 
theory at low energies. 

The method usually used to calculate perturbative QCD contributions to quantities such 
as R T consists in rewriting the original expression, which involves integration over small 
values of momentum, by means of Cauchy's theorem and thereby convert it into a contour 
integral in a region of relatively large momentum. But the occurrence of incorrect analytic 
properties of the perturbative approximation, such as the appearance of a ghost pole and 
an unphysical cut, make it impossible to exploit Cauchy's theorem in this manner. 

In this paper we will use analytic perturbation theory which provides a correlator 
II(g 2 ) with the correct analytic properties. For example, the one-loop analytic running 
coupling constant a an (g 2 ) in this approach is given byQ 



where a = a s j-K and 6/2 is the first /3-function coefficient. It is clear that, unlike the standard 
perturbative expression, this one has no unphysical ghost pole and, therefore, possesses the 
correct analytic properties, arising from the Kallen-Lehmann analyticity that reflects the 
principle of causality. The nonperturbative term which appears in the running coupling 
constant does not change the ultraviolet limit of theory. A distinguishing feature of this 
approach is the presence of a universal limit at q 2 = 0, independent of both the scale 
parameter A and choice of renormalization scheme. This quantity turns out to be stable 
with respect to higher loop corrections. 

In order to compare theoretical results with experimental measurements one has to esti- 
mate the uncertainty connected to the renormalization scheme dependence of the perturba- 
tive approximation. It is well known that at low energies this ambiguity becomes consider- 
able (see, e.g., Ref. Analytic perturbation theory improves the situation and gives very 
stable results over a wide range of renormalization schemes, as has been demonstrated for 
the e + e~ annihilation ratio in Ref. |4]]. 

We will perform the calculation of R T at the two- and three-loop levels and study the 
renormalization scheme dependence within analytic perturbation theory. All results will 
be compared with standard perturbative ones. This investigation is a continuation of the 
program initiated in Ref. 0. 



We use the definition q 2 < in the Euclidean region. 




2 T 1 A 2 



(1) 



b [ln(-g 2 /A 2 ) A 2 + q 2 
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II. PROBLEM OF QCD PARAMETRIZATION OF R T 



The ratio of hadronic to leptonic r-decay widths 

_ T(t~ -> z/ T hadrons) 
T(r -> i/ T < v{) 

has the following theoretical representation 

R T = 3 S EW (IK/ + |KsD(l + 5qcd), (3) 

where V u d and V^ s are the CKM matrix elements, Sew is the electroweak factor, and £qcd 
is the pure QCD correction [H], which we will calculate perturbatively. It is convenient to 
introduce some definitions: Im II ~ 1 + r for the perturbative part of correlator II (q 2 ) and 
D ~ 1 + d for the perturbative part of Adler function D(q 2 ). Then one can write <5qcd as 
an integral over timelike momentum s: 

If one uses the standard perturbative approximation this integral is ill-defined because of 
unphysical singularities within the range of integration; standard perturbation theory (PT) 
fails. The most useful trick to rescue the situation is to appeal to analytic properties of the 
hadronic correlator II(g 2 ) ||. This opens up the possibility of exploiting Cauchy's theorem 
by rewriting the integral in the form of a contour integral in the complex g 2 -plane with the 
contour being a circle of radius M 2 : 

However, as noted, the transition to the contour representation @ requires certain analytic 
properties of the correlator, namely, that it must be an analytic function in the complex 
g 2 -plane with a cut along the positive real axis. The correlator parametrized by the running 
coupling constant in standard perturbation theory does not have this virtue J7||| . Moreover, 
a renormalization group analysis gives a running coupling constant determined in the space- 
like region, whereas the initial expression (|J) for R T contains an integration over timelike 
momentum. Thus, we are in need of some method of continuing the running coupling con- 
stant from the spacelike to the timelike region that takes into account the proper analytic 
properties of the running coupling constant. Because of this failure of analyticity Eqs. ([|) 
and (|J) are not equivalent in PT. 

Nevertheless, if one uses Eq. (|5|) with the PT approximation for d(z), thus calculating 
the contour integral in a naive way, a new problem will emerge. It is known that the lower 
the energy scale we consider, the greater the uncertainty there is related to the choice of 
renormalization scheme (RS). If one turns to the _D(g 2 )-function, which can be extracted 
from the ratio R e + e - for e + e~ annihilation, one finds that the RS dependence of theoretical 
predictions is considerable 0]. At the energy scale around M T this uncertainty increases 
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drastically. The transition to the contour integral softens this ambiguity, but it remains 
large. Moreover, as we noted, within PT this transition to the contour integral representation 
cannot be performed self-consistently. Therefore, there is no consistent method of calculating 
the inclusive decay of the r-lepton into hadrons in the framework of standard perturbation 
theory. 

Recently, a new approach that eliminates the two problems described above has been 
proposed This is the so-called analytic perturbation theory (APT), in which the 

hadronic correlator II(g 2 ) has the correct analytic properties. In this approach, the <i an (g 2 ) 
function, which is analytic in the cut g 2 -plane, can be expressed in terms of the spectral 
density p(cr), the basic quantity in APT (to distinguish APT and PT we have introduced 
subscripts "an" and "pt" instead of QCD) 

1 r°° dn 

4n(g 2 ) = -/ -^sfa). (6) 

it Jo a — q z 

We can write down some useful dispersion relations, which are consequences of the ana- 
lytic properties of the correlator. The first one is 

Jo [s — q 2 ) 
which is inverted by the following formula, 




Here, the contour lies in the region of analyticity of the D(g 2 )-function. Then let us rewrite 
the last expression in terms of the spectral density p(a) f|: 

. . 1 f 00 da , . 
r**(s) = - / —p(a). (9) 

71 Js <J 

If we substitute the last formula into integral (|j), we can easily verify that the contour 
representation @ holds [|]. Using the spectral representation flU), we can rewrite 5 an in 
terms of p(cr): 




The first term in this formula is simply the value d an (0). Once an expression for p(a) is 
given, we can calculate the integral (|10|). 

Consider <5 an at the three-loop level. The running coupling constant a = a s /n satisfies 
the NNLO renormalization group equation: 

P 2 -^ = ~a 2 (^ + cia + c 2 a 2 ). (11) 
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In the modified minimal subtraction scheme, MS, for three active flavors, we have b = 
9/2, d = 16/9 and cf s = 3863/864. The renormalization- group improved perturbation 
expansion for d pt (q 2 ) is given byf\ 

d P t(q 2 ) = a(q 2 )[l + d ia (q 2 ) + d 2 a 2 (q 2 )]. (12) 

In the following, values of d\ and d 2 will be needed; in the MS-scheme they are d^ = 1.6398 
and df* = 6.3710 f|. 

In the APT approach, the spectral density is defined as the imaginary part of the per- 
turbative approximation to d pt {q 2 ) on the physical cut: 

p{<?) = Qo(o) + dtQi (a) + d 2 g 2 (cr), (13) 

where 

g n (a) = lm[a n+1 (a + te)]. (14) 

The first of these spectral densities gives an analytic expression for the running coupling 
constant a an (g 2 ): 




This is what was called the spacelike running coupling constant in Ref. ||. In the one- loop 
approximation it leads to Eq. (]!]). 

For the following considerations, it is important that d^q 2 ) and a an (g 2 ) have a universal 
limit at the point q 2 = 0. This limiting value, generally, is independent of both the scale 
parameter A and the order of the loop expansion being considered. Because d an (0) and 
Oan(0) are equal to the reciprocal of the first coefficient of the QCD /^-function, they are also 
RS invariant (we consider only gauge- and mass- independent RSs). The existence of this 
fixed point will play a decisive role in the very weak RS dependence of our results. 

The next point we should like to note is that the two approaches, APT and PT, coincide 
with each other in the asymptotic region of high energies. The procedure of implementing 
analyticity described above does not change the ultraviolet limit of the theory, i.e., d^q 2 ) 
has the same asymptotic behavior as q 2 — > — oo as the standard PT Adler function. 

Therefore, the above gives a self-consistent framework in which we can estimate the 
-R r -ratio and extract from the measured value of R T the QCD scale parameter A. 



2 We have made a few changes in notation from that given in Ref. |5): Now a = a s /7T, and 
consequently we have denoted by d\ and d 2 what we called d 2 and d% previously, apart from a 
power of 4. 
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III. RENORMALIZATION SCHEME DEPENDENCE 



In this section we will investigate the RS dependence of the APT approach.^ The Adler 
function d an , as in PT, is parametrized by a set of RS parameters, but only two of them 




are independent. There is a RS invariant combination u 2 |yj °f parameters that binds 
them: 

oo 2 = c 2 + d 2 - c x d x - dx, (16) 

which in our case equals 5.2378. Here, C\ is RS invariant in the class of schemes mentioned 
above. Thus, we can choose d\ and c 2 as independent variables, which define the RS, and 
parametrize the QCD correction in terms of them. 

To find G? an (g 2 ) we use Eqs. (pif), flI3|), and ([]) and solve the transcendental equation for 
the running coupling constant on the physical cut lying along the positive real axis: 

-di + i + Cilnf^-) + F«(a), (17) 
where at the two-loop order 

and at the three-loop order 

^M^'M + of {l + cix){l d ; cix + c2x2y (19) 

These expressions hold in any renormalization scheme. 

It is necessary to extract from the fit in the MS-scheme and after this procedure 
to pass on to other renormalization schemes with the RS-invariant parameter A^jg. To do 
the fit we took the values V ud = 0.9752, V us = 0.2205, S EW = 1-0194 and M T = 1.777 GeV 
T2],|T3"||. For numerical estimations, we used the world average value R T = 3.633 ± 0.031 



14], which leads to the following value of the QCD scale parameter in three-loops: A^ 



871 ± 155 MeV. The corresponding value of 4n(M 2 ) is 0.1446 ± 0.0088 and of a an (M 2 ) is 
0.3962 ±0.0298. To compare this result with the PT prediction we did the same fit using the 
contour integral representation (§) and got A^ = 385±27 MeV, d pt (M 2 ) = 0.1339±0.0069, 
and a p t(M 2 ) = 0.3371 ± 0.0141. Thus, the values of the Z)-function are close to each other, 
but the value of A^ in the APT approach is much larger than the one obtained in PT. As 
in Ref. ||, we found that this value is very sensitive to the value of R T . In other words, <5 an 
depends so slightly on that a 0.9% error in R T gives 18% error in the value of A^g. We 
illustrate this feature in Table 1. According to the table, when we change M 2 /A 2 from 2.0 
to 6.5 (corresponding to a variation of A from 1.256 GeV to 0.697 GeV), 5 an is only altered 
by about 20%. The sensitivity to Aj^g increases as M 2 /A 2 gets smaller. 



3 In the framework of the conventional approach, this problem has been studied in Ref. [10|. 



6 



For comparison, we also include the two-loop result: A|2g = 918 ± 151 MeV, <i an (M^) = 
0.1431 ± 0.0084, and the running coupling constant a an (M T 2 ) = 0.3983 ± 0.0236. The 
perturbative result is = 458 ± 31 MeV, d pt (M T 2 ) = 0.140 ± 0.0067, and a pt (M 2 ) = 
0.3544 ±0.0157. 

The reader should further note that the convergence properties of the APT series are 
much better than those of the PT case; in the improved scheme, the main contribution comes 
from the diagrams of lowest order. For example, the value of the three-loop APT correction 
to d an (M^) is smaller than the corresponding PT contribution to d pt (M^) by about a factor 
of four. 

Having found A^jg, we will now study how 5 an varies with a change of renormalization 
scheme. A natural way to do so is to supplement results in a certain scheme with an estimate 
of variability of predictions over a whole range of a priori acceptable schemes. To restrict 
this uncertainty one has to impose some requirement eliminating schemes with unnaturally 
large coefficients that would introduce large cancellations. It was proposed to calculate 
theoretical predictions over the set of schemes satisfying the condition Jl5| 



|c 2 | + \d 2 \ +cx\di\ + di <2\u 2 \, (20) 

which allows schemes with a degree of cancellation the same or smaller than what occurs in 
the scheme obeying the so-called principle of minimal sensitivity (PMS). 

In the MS-scheme we adopt <5 ariiPt = 0.1881. Let us consider some points belonging to 
the domain described above and compare corresponding values of 5 pt . For example, suppose 
we choose the points A = (—1.6183,0) and B = (0.9575,0), where the first coordinate is 
d\ and the second is c 2 . Both points lie on the boundary of the domain, i.e., they have 
the same degree of cancellation as does the PMS scheme. But S pt (A) = 0.1999, while 
Spt(B) = 0.1886. Therefore, even with such a restriction imposed on RSs we have about a 
5% deviation of PT results from each other. The difference between APT results is much 
smaller: 5 an (A) = 0.1866 and S an (B) = 0.1881, so we have only 0.8% deviation. For clarity 
we display our three-loop results in the form of a contour plot, in Fig. |I[ Two-loop results 
are 5 pt (A) = 0.2182 and S pt (B) = 0.2070. So we have again about a 5% variation in the PT 
results. In APT one has S ari (A) = 0.1899 and 5 an (B) = 0.1890, which corresponds to only a 
0.4% variation. 

Note that the MS-scheme itself does not belong the domain (|20"D . Therefore, it is worth- 
while to consider some schemes lying outside the domain. In Ref. |10[ it was shown that the 



so-called V^-scheme [^6| lies very far from the domain described above and gives so a large 
value of 5 pt that it cannot be used at this low energy. For the ^/-scheme we have d\ = —0.109 
and C2 = 26.200. The three-loop PT result is 5 pt = 0.3115, corresponding to about a 66% 
deviation from the MS-scheme. On the other hand, if we turn to APT we have S an = 0.1878, 
i.e., only about a 0.2% deviation from the MS-scheme. So the ^-scheme is still useful at 
this energy in APT. 

This stability of the APT method is a consequence of the existence of the RS-invariant 
fixed point, (i an (0). In PT at high energies the weak RS dependence is a consequence of the 
small value of the coupling constant. At lower energies uncertainty increases. In APT, at 
high energies, the situation is the same, but at low energies the theory has a universal limit, 
which restricts the RS ambiguity over a very wide range of momentum. 
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Another way to illustrate the remarkable stability of APT is to calculate the spectral 
functions g n (cr) given by Eq. fll4]) ; one sees that Qi{a) is much smaller than Qo((t) over the 
whole spectral region. The same statement is true for the relationship between Qi(cr) and 
^ 2 (o")- This monotonically decreasing behavior reduces the RS dependence strongly, since 
the perturbative coefficients d\ and d>2 in expression ( |I~3D for p(a) are multiplied by these 
functions. 

IV. CONCLUSION 

We have considered inclusive r-decay in three-loop order within analytic perturbation 
theory, which, in contrast to standard perturbation theory, is a self-consistent procedure. 
We summarize the following important features of this method. 

First, the method maintains the correct analytic properties and leads to a self-consistent 
definition of the procedure of analytic continuation. 

Second, it gives an unusually large value of the QCD scale parameter A, connected to 
the presence of nonperturbative contributions that appear in the APT method. The value 
of A is very sensitive to the precise value of R T . 

Third, the three-loop correction is smaller than the one found in standard perturbation 
theory. Hence, the APT approach is stable with respect to higher loop corrections. 

Fourth, the RS dependence of the results obtained is reduced drastically. For example, 
the V^-scheme, which gives a very large discrepancy in standard perturbation theory, can be 
used in analytic perturbation theory without any difficulty. Thus, the APT predictions are 
practically RS independent over a wide region of RS parameters. 

Therefore, we have a self-consistent method within analytic perturbation theory of cal- 
culating semileptonic r-decay. 

In this paper, we have not considered the standard nonperturbative power corrections. 
The process of enforcing analyticity changes the perturbative contributions by incorporating 
some nonperturbative terms; consequently, power corrections should also be changed. The 
role of power corrections is still not understood within the APT approach, and requires a 
separate discussion. We hope to accomplish this in our subsequent papers. 

ACKNOWLEDGEMENTS 

The authors would like to thank D.V. Shirkov and O.P. Solovtsova for interest in this 
work and for useful comments. Partial support of the work by the US National Science 
Foundation, grant PHY-9600421, and by the US Department of Energy, grant DE-FG-02- 
95ER40923, is gratefully acknowledged. The work of ILS is also supported in part by the 
University of Oklahoma, through its College of Arts and Science, International Programs, 
Vice President for Research, and the Department of Physics. 



8 



REFERENCES 



D. V. Shirkov and I. L. Solovtsov, JINR Rep. Comm. 2[76]-96 5 (1996) , |hep-ph/ 9604363 

D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79, 1209 (1997). 
P. A. Raczka and A. Szymacha, Phys. Rev. D 54, 3073 (1996). 

I. L. Solovtsov and D. V. Shirkov, |hep-ph/9711^BT| . 

K. A. Milton, I. L. Solovtsov and O. P. Solovtsova, Phys. Lett. B 415, 104 (1997). 

E. Braaten, Phys. Rev. Lett. 60, 1606 (1988); 

E. Braaten, S. Narison and A. Pich, Nucl. Phys. B 373, 581 (1992). 

H. F. Jones and I. L. Solovtsov, Phys. Lett. B 349, 519 (1995); 

H. F. Jones, I. L. Solovtsov and O. P. Solovtsova, Phys. Lett. B 357, 441 (1995). 

K. A. Milton and I. L. Solovtsov, Phys. Rev. D 55, 5295 (1997). 

S. G. Gorishny, A. L Kataev and S. A Larin, Phys. Lett. B 259, 144 (1991). 

P. A. Raczka, |hep-ph/ 9707366 . 



P. M. Stevenson, Phys. Lett. B 100, 61 (1981); 
Phys. Rev. D 23, 2916 (1981). 

E. Braaten and Chong Sheng Li, Phys. Rev. D 42, 3888 (1990). 
W. J. Marciano and Sirlin, Phys. Rev. Lett. 61, 1815 (1988). 
Particle Data Group, Phys. Rev. D 54, 77 (1996). 
P. A. Raczka, Z. Phys. C 65, 481 (1995). 
W. Fischler, Nucl. Phys. B 129, 157 (1977); 

T. Appelquist, M. Dine and I. J. Muzinich, Phys. Lett. B 69, 231 (1977). 



9 



FIGURES 
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TABLES 



M r 2 /A 2 e NL ° m t va 2 e° <e L ° 

2.0 0.2090 0.2106 4.5 0.1820 0.1857 

2.5 0.2016 0.2039 5.0 0.1785 0.1824 

3.0 0.1955 0.1983 5.5 0.1753 0.1795 

3.5 0.1904 0.1935 6.0 0.1724 0.1767 

4.0 0.1859 0.1894 6.5 0.1698 0.1743 



TABLE I. NLO and NNLO predictions for 5 an in the MS scheme. 
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